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Abstract 

Asymptotic behavior of the anomalous dimensions of Wilson operators with high spin and 
twist is governed in planar J\f = 4 SYM theory by the scaling function which coincides at strong 
coupling with the energy density of a two-dimensional bosonic 0(6) sigma model. We calculate 
this function by combining the two- loop correction to the energy density for the 0(n) model with 
two-loop correction to the mass gap determined by the all-loop Bethe ansatz in A/" = 4 SYM 
theory. The result is in agreement with the prediction coming from the thermodynamical limit of 
the quantum string Bethe ansatz equations, but disagrees with the two-loop stringy corrections 
to the folded spinning string solution. 



1 Introduction 



The AdS/CFT [Tj establishes the correspondence between Wilson operators in TV = 4 super 
Yang-Mills theory (SYM) and states of strings spinning on AdSs x background [21 [3]. It has 
been recently recognized that there exists a remarkable relation between both theories and the 
two-dimensional bosonic 0(6) sigma model. This relation emerges when one studies anomalous 
dimensions of Wilson operators in planar A/" = 4 SYM theory at strong coupling in the limit [1] 
when the Lorentz spin of the operators grows exponentially with their twist L 

3 = = fixed , N, L^oo. (1.1) 

InA^ ^ ^ 

In this limit, the anomalous dimensions grow logarithmically with A^ and their leading asymptotic 
behavior is governed by the scaling functioE0 depending on j and the t' Hooft coupling = 
gyj^Nc/ (4:71)'^ . Using the dual description of Wilson operators as folded strings spinning on 
AdSs X and taking into account the one-loop stringy corrections to these states [5j, Alday 
and Maldacena [6] conjectured that the scaling function should coincide at strong coupling with 
the energy density of a two-dimensional bosonic 0(6) sigma model. Recently, this relation was 
established in planar A/" = 4 SYM theory at strong coupling \7\ using the conjectured integrability 
of the dilatation operator |8|. 

The 0(6) sigma- model and, in general, two-dimensional bosonic 0(n) sigma-models are 
among the best studied field theory models. Their popularity is explained by the fact that they 
can be used as low dimensional toy models of QCD: they are asymptotically free in perturbation 
theory, their classical conformal invariance is broken and the mass of the physical excitations is 
dynamically generated by the dimensional transmutation mechanism [9l [10]. In addition they 
are integrable and many physical quantities are exactly calculable [11]. These models can also be 
studied non-perturbatively with a help of lattice Monte Carlo simulations making use of a very 
efficient simulation algorithm (the cluster algorithm). 

Using the known exact S-matrix of the 0(n) sigma-model, a linear thermodynamical Bethe 
ansatz (TBA) integral equation can be derived describing the free energy in the presence of an 
external field coupled to one of the Noether currents of the modelH The free energy can also 
be calculated in perturbation theory for large values of the external field due to asymptotic 
freedom. The original motivation of this calculation [12] was that by comparing the two results 
for the free energy one can calculate the physical mass m of the 0(n) particles in terms of the 
perturbative A-parameter (dynamical scale defining solutions to Gell-Mann-Low equation). The 
exact mass gap m/A determined this way has been checked by using Monte Carlo results and 
also by comparing finite volume mass gap values computed in perturbation theory [13] to those 
obtained from a (nonlinear) TBA integral equation 1^^. 

In this paper, we employ tools developed for the 0(n) sigma-model to compute the scaling 
function in the AdS/CFT. 

1.1 Scaling function in = 4 SYM 

The Wilson operators under consideration are built from L complex scalar fields and A^ light-cone 
components of the covariant derivatives. Their minimal anomalous dimension has the following 

"'^More precisely, for given spin N and twist L, anomafous dimensions of Wilson operators occupy a band. The 
scaling function describes the minimal anomalous dimension in this band. 

^The energy density is related to the free energy through Legendre transformation. 
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behavior both at weak and at strong couphng [U El [8], [5], [15] 



-iNAa) = [2r,,sp(^7) + e(^7,j)]lniV + ... , (1.2) 

where elhpses denote terms suppressed by powers of 1/L. Here the first term inside the square 
brackets has a universal, j— independent form and it involves the cusp anomalous dimension [161 
[T7] . The dependence on the twist resides in e{g,j) which is a nontrivial function of the 't Hooft 
coupling and the scaling variable j normalized as e{g, 0) = 0. 

At weak coupling, the scaling functions Tcusplfi') and e{g,j) can be found in a generic (su- 
persymmetric) Yang-Mills theory in the planar limit by making use of the remarkable property 
of integrability [H [IT]. In maximally supersymmetric Af = 4 theory, these functions can be 
determined in the planar limit for arbitrary values of the scaling parameter j and the coupling 
g as solutions to BES/FRS equations proposed in El [H]. These equations predict the scaling 
function to be a bi-analytical function of g^ and j.u At strong coupling, the asymptotic ex- 
pansion of the cusp anomalous dimension in powers of 1/g was derived in [20l [211 [221 [231 [24] . 
It turned out that this expansion suffers from Borel singularities [21] indicating that rcusp(fi') 
receives nonperturbative corrections at strong coupling defined by the scale m ~ g^^^e''^^ j6l[2^. 

At strong coupling, the scaling function €{g,j) has a more complicated form and its properties 
depend on a hierarchy between g and j. As was shown in [7], for g ^ oo and j/m = fixed (with 
m given by (11 ■4p below), the FRS equation for the scahng function e{g,j) can be casted into a 
form identical to the TBA equations for the nonlinear 0(6) sigma model [12]. This leads to the 
identification of the scaling function as the energy density £0(6) in the ground state of the 0(6) 
model corresponding to the particle density p 

£o(6) = ^ , P^2' 

As was already mentioned, the 0(6) sigma model has a nontrivial dynamics in the infrared and 
it develops a mass gap. Remarkably enough, similar phenomenon also occurs for the scaling 
function in A/" = 4 theory at strong coupling. Namely, the scaling function e{g,j) depends on a 
new dynamical scale [3, [25] 

23/47rV4 

m = kg'/'e'-^[l + Oil/g)], ^ = ^^'^^ 

and the same scale fll.4p defines nonperturbative corrections to the cusp anomalous dimension at 
strong coupling [211 [26]. Later in the paper, we will compute subleading corrections to the mass 
scale (11. 4p . The relation (II. 3p holds for g ^ 00 with j/m = fixed and the scale (11.40 is identified 
as the mass gap of the 0(6) model. 

1.2 Scaling function in AdS/CFT 

The AdS/CFT correspondence relates the anomalous dimension (11.21) at strong coupling to the 
energy of a folded string spinning on the AdSs x background [21 [3] 

A = N + L + ^nA9), (1-5) 



^The properties of this function at weak coupling and large j were studied in [19j 
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with N and L being angular momenta on AdSa and S^, respectively. Semiclassical quantization 
of this stringy state yields the expansion of the anomalous dimension 'jN,L{g) in powers of 1/g. 
In agreement with fll.2p . the coefficients of the expansion scale logarithmically with in the 
limit (11. ip and determine stringy corrections to the scaling function. 

For the cusp anomalous dimension T cuspid), the first three coefficients of 1/g expansion were 
computed in Refs. [21 |5l ^7\. As a nontrivial test of the AdS/CFT correspondence, they were 
found to be the same as in the A/" = 4 SYM theory plj0 ' For the scaling function similar 
calculation was performed by two different approaches - from two-loop stringy corrections to the 
folded spinning string solution [28] and from thermodynamical limit of quantum string Bethe 
ansatz equations [151 ISH [30j, leading to (in notations of [28] ) 



eig,j)+j = 2f 



g + 



TT 



002 

2 



31n£ + 4(ln£)2) + O (l//) 



+ Oit), (1.6) 



with i = j / (4(7). This relation is valid at strong coupling for j <^ g, or equivalently £ <^ 1. Here 
the first two terms inside the square brackets describe, correspondingly, the classical expression 
and one-loop correction to the scaling function e{g,j). The last term describes the two-loop 
correction and it depends on a constant go2- The two approaches mentioned above produce two 
different values of go2 



qo2 



re£p5) 



-2K - -ln2 + -, 
2 4 ' 



qo2 



refgn] 
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2 



(1.7) 



with K being the Catalan constant. The two results agree with each other in term ~ ln2 but 
disagree in the rest. The reason for this discrepancy remains unclear. 

The semiclassical approach allows us to calculate the scaling function (11. 6p in the form of a 
double series in 1/g and i"^. It does not take however into account nonperturbative corrections 
to the scaling function which are exponentially small as g ^ oo. Alday and Maldacena ^ put 
forward the proposal that the scaling function e{g, j) can be found exactly at strong couphng in the 
limit j g and j/m = fixed (with m defined in (II. 4p ). They argued that quantum corrections in 
the AdSs x sigma model are dominated in this limit by the contribution of massless excitations 
on whose dynamics is described by a (noncritical) two-dimensional bosonic 0(6) sigma-model 
equipped with a UV cut-off determined by the mass of massive excitations. In terms of parameters 
of the underlying AdSs x sigma model, the exact value of the mass gap is given by (II. 4p . The 
dependence of the mass scale (II. 4p on the coupling constant is fixed by the two-loop beta-function 
of the 0(6) model whereas the prefactor k was determined in [6j by matching the first two terms 
of the semiclassical expansion (ll.6p into known one-loop perturbative correction to the energy 
density of the 0(6) model. 

The subleading 0(1/ g) correction to the scaling function (II. 6p involves both constant term 
go2 and logarithmically enhanced terms. As was shown in [6], the latter terms are controlled 
by renormalization group. It is the constant go2 that lies at crux of the relation (II. 6p to two- 
loop order. The relation (II. 3p allows us to determine this constant by computing corrections to 
the energy density of the 0(6) model and by matching the resulting expression for the scaling 
function into (11.60 . 



Note that the semiclassical approach does not take into account nonperturbative corrections to rcusp(<?) 
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1.3 Two-dimensional 0(6) sigma model 

Substitution of fll.6p into relation (11 .Sp yields a definite prediction for the energy density of the 
0(6) model in the regime of large particle density p ^ m. In this regime the model is known to 
be weakly coupled and the energy density can be computed perturbatively. 

The two-dimensional 0(6) sigma model is an exactly solvable quantum field theory. It is 
asymptotically free at short distances, while in the infrared it develops a mass gap. The massive 
excitations form the vector multiplet of the 0(6) group and their S-matrix has a factorized form 
|llj . This makes it possible to calculate the mass gap m in terms of the parameter A. The idea is 
to consider the 0(6) model in the presence of an external field h coupled to one of the conserved 
charges, say Q^"^ (see Eq. (13. ip below), and calculate the change in the free energy density in two 
different ways: from thermodynamical Bethe ansatz and from perturbative expansion [12] . 

If the external field exceeds the mass gap, h > m, a. finite density of particles p is formed in 
the ground state. The corresponding ground-state energy density £o(6)(p) is given by 

/B pB 
— x(^)cosh^, p = / —x{0), (1.8) 

where the rapidity distribution satisfies an integral TBA equation 

j-B 

X{0)= d6'K{9-9')x{9')+m cosh 6. (1.9) 

J-B 

Here the kernel K{9) = (log S'(6') )'/(27ri) is given by the logarithmic derivative of the scattering 
matrix of the particles corresponding to the largest eigenvalue of the charge Q^^ 



K{9) 



271 J V 27ry V2 27ry V2 27ry cosh^ 



;i.io) 



where 'il){x) = r'(x)/r(x) denotes the logarithmic derivative of the gamma function. 

The free energy density J^{h) can be obtained from (I1.8P through the Legendre transformation, 
J^{h) = miup [e{p) — hp]. Due to asymptotic freedom, for h ^ m the change in the free energy 
density can also be calculated perturbatively with the result 

- m = --l--l-^ + 0{a')] , (1.11) 
where the coupling a = a{h) is defined in a renormalization group invariant way by 

- + -lna = ln- , 1.12 

with Aj;jg being the A scale in the MS scheme. The first two terms in the right-hand side of (II. lip 
were found in [12j while calculation of the coefficient in front of a is one of the main results of 
this paper. 

At large h, the solution to the TBA equations (11. 8p and (11. 9p can be constructed using 
the generalized Wiener- Hopf technique [12j. The result for J-'{h) — J-'{0) is given by series in 
[ln(/i/m)]~^. The calculation is rather involved and only the first few terms have been calculated 
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so far |12j. Matching these terms into perturbative expansion (ll.lip . one can estabhsh the 
relation [12] between the mass gap m and the scale Aj^ 



C = ln 



-(31n2-l)-lnr(f). 



(1.13) 



■MS 



Taking into account this relation, we find from (11.111) and (I1.12p the following expression for the 
ground state energy density as a function of p/m 



The relation (11.141) defines the energy density eo{6) as a function of the particle density p and the 
mass gap m. We can now apply (11.31) to translate it into the dependence of the scaling function 
^{QiJ) oil 3 and m in the limit m <^ j <^ g. 

To obtain the strong coupling expansion of e{g,j) from (II. 3p we also need the explicit form 
of the function m = m{g). To leading order it is given by (ll.4p while the subleading correction 
to m will be computed below (see Eq. (I2.4p ). Replacing the mass scale m in (11.140 by its explicit 
expression (11.40 . we find that the scaling function admits the same perturbative expansion as 
(I1.6P and provides a definite prediction for the constant go2- This constant depends on subleading 
corrections both to the energy density £0(6) and to the mass scale m. We find that with these 
corrections taken into account the relation (11.30 leads to go2 = — 31n2/2 + 11/4 (see Eq. (II. 7p ). 
in agreement with the quantum string Bethe ansatz result of [30] . 

The paper is organized as follows. In Sect. 2 we summarize the properties of the scaling 
function in planar J\f = 4 SYM at strong coupling. We show that for m j <^ g the scaling 
function has the same form as (II. 6p and obtain the expression for the constant go2 which involves 
corrections both to the mass gap m and to the energy density £0(6)- Then, we compute subleading 
0(1/(7) correction to the mass scale (II. 4p . In Sect. 3, we calculate the ground-state energy density 
of the 0(n) sigma model for large particle density p using perturbation theory technique. We 
obtain an expression for the energy density 6o{n) to second order in the effective coupling constant 
expansion and demonstrate that it is in agreement with numerical solution of the TBA equations. 
Finally, we use the obtained expressions for m and £0(6) to compute the constant go2- Section 4 
contains our concluding remarks. Some technical details are presented in two Appendices. 

2 Scaling function in A/^ = 4 SYM at strong coupling 

For g —>■ oo and j/m = fixed, the scaling function e{g,j) is related to the energy density of the 
0(6) model through relation (11.30 . To compute e{g,j) from (11.30 we have to accomplish two 
tasks: (i) to determine the energy density eo(6) as a function of particle density p = j/2 and 
mass gap m, and (ii) calculate the mass scale m from the FRS equation as a function of the 
coupling constant g. 




(1.14) 



where a = a{p) denotes another useful coupling 



— lna = ln hln — hC- 

a 4 m 2 



(1.15) 
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The leading order expression for the mass scale is given by (]1.4p . The subleading correction 
to m can be parameterized as follows 



m 



kg 



1/4 Q-ng 



ng 



(2.1) 



with k defined in (11.41) and (/—independent parameter mi to be determined. For the energy 
density £o(6)) we make use of the results of previous studies of the TBA equations for the two- 
dimensional 0{n) model [12]. For arbitrary j/m, solutions to the TBA equations (11.81) and 
(11.91) do not admit a simple analytical representation. However, they can be constructed in the 
limits j/m <^ 1 and j/m 3> 1, which correspond, respectively, to (nonperturbative) small and 
(perturbative) large particle density regimes: 

• For j <C m, the integral equations (II. 8p and (11.91) can be solved by iterations leading 
to [HI El [71 



m 24 



+ . 



(2.2) 



where expansion runs in powers of j/m. This regime was extensively studied both nu- 
merically and analytically and subleading corrections to (12.21) were recently computed in 



For j ^ m, the expression for the scaling function follows from the known perturbative 
result for the energy density of 0(6) model [121 El [7] 



where In k 



\ ln(j'/m) 



3ln(Kln(j/m)) + i 9 ln^(Kln(jym)) + ei 



In (j'/m) 



+ 



16 



ln^(j/m) 



+ . . . 



(2.3) 



I — I ln2 — I lnr(|) and the constant ei remains unknown. 



The relation (12. 3p resums through renormalization group (an infinite number of) perturbative 
corrections m. 1/ g which are proportional to and are enhanced by powers of \n{j/g) [HI EB]- 
Replacing the mass scale in (12. 3p by its expression (12.11) and re-expanding the right-hand side of 
(12. 3p in powers of 1/ g, we arrive at the relation (I1.6P and obtain the constant go2 as 



qo2 



9 9 
- + 8mi + -ei. 



(2.4) 



It depends on the parameters mi and ei entering (12. ip and (12. 3p . respectively. In this section, 
we compute mi and return to ei in section 3. 



2.1 Mass scale 



The dependence of the mass scale m on the coupling follows univocally from the FRS equation. 
As was shown in [7], it has the following form 



m 



8^2 



-^e-^Re 

TT 



dt e^(*-^/4) 
t + mg 



(r+(t) + ^r_(t)) 



(2.5) 



where T^{t) are real functions of t also depending on the coupling constant but independent on 
the scaling variable j. In addition, these functions have a definite parity, r_|_(— t) = ±r_|_(t), and 



6 



satisfy an integral equation [24] which follows from the BES equation [18] . To save space, we do 
not present it here and refer an interested reader to [3, [21] for details. 

It is interesting to note that the functions T^(t) also play a distinguished role in determination 
of the cusp anomalous dimension in planar Af = 4 SYM theory. Namely, for arbitrary coupling 
rcusp(5') can be derived from the asymptotic behavior of T_^{t) at small t 



TcusM = -^9 (r+(o) + zr_(o)) 



(2.6) 



This relation was used in [21] to obtain the asymptotic expansion of Fcusplfi') in powers of 1/g. 

The functions T±{t) were constructed in in the form of Neumann series over Bessel 
functions 



r_(t) 



fc=0 

oo 

Hi 

k=0 



\k+l 



\k+l 



J2k{t)T2k{g) , 



J2k-lit)T2k-l 



(2.7) 



Here the dependent expansion coefficients are given by 

r-1 = 1 , r,(,) = -Irf ) + 1 [c;(,)rf + 4(,)rf 



{2p) 



p=l 



(fc>0). 



f2.8) 



where the dependence on k is carried by the coefficients defined as 



T{k + p 



2k-l 



(2.9) 



r(fc + i)r(i)' ' ' r(A; + i-p)r(2; 

while the dependence on the coupling resides in Cp{g). 

The functions T±(t) defined in (12. 7p are uniquely specified by the expansion coefficients Cp{g). 
The latter satisfy the quantization conditions 



J2{2nsrc;ig) 



Tip 



p>0 



2T( 



r(|)r(i-.) 
r(!-^) 



0(1/^?) 



(2.10) 



J^i^nsr [c,{g) + 4{g){2p - |) - 1)] 



r(i)r(i 



p>0 



r(i-^) 



+ 0{l/g) 



with Cn 



'2' '^0 



and s being arbitrary. Comparing coefficients in front of powers of s in 
both sides of these relations, we can determine c^{g) in the form of asymptotic series in 1/g. In 
this manner, we get 



ctig) = -3 — 

TT 

, 3 ln2 

Cl (^) = 7 — - 
4 77 



1 

4 



+ 0{l/g) 
^0{l/g), 



(2.11) 



Substituting these relations into (12. 7p and (12. Sp and performing summation over k in the right- 
hand side of (12. sp we obtain after some algebra 



TMt)+iT_iit) = -Voit)-iA7rg) 



-1 



|-31n2) tVoit) 



31n2 



Viit) 



+ O (1//) , (2.12) 
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where the notation was introduced for the functions 

1/4 



TT J_i \1 — U J 1 + U 

The first term in the right-hand side of (12.121) coincides with the leading-order solution found in 
j21j . Moreover, taking t = in fl2.12p we verify using (12.61) that 

reusp(^)/(2<7) = -r+(o) - ^r_(o) = 1 - ^ + . . . , (2.i4) 

in agreement with the known strong coupling expansion of the cusp anomalous dimension. 



2.2 Correction to the mass gap 



Let us now apply (12. 5p and (12. 7p to compute the mass scale at strong coupling. We observe 
that the t— integral in the right-hand side of (12.50 receives a dominant contribution from the 
region t ^ g. Trying to apply (I2.12p we recognize that the expansion (I2.12p is not well-defined 
in this region because, due to the presence of t inside square brackets, expansion parameter is 
t/{4:7rg) = 0{g^). Thus, in order to compute the mass scale from (12. 5p . we have to resum the 
whole series (I2.12p in the double scaling limit t,g oo with t/g = fixed. Fortunately, this 
particular limit was already studied in [21]. 

In the double scaling limit, it is convenient to change the integration variable in (12. 5p as 
t — *• Airgit and expand the function r^{4n git) into series in 1/g with t = 0{g^). We find from 
(I2.12p that the expansion has the following form 



T^Angzt) + iT_{47rgit) = fo{t)Vo{4n gt) + f,{t)Vi{4ngt) 



(2.15) 



with /o(t) = -1- (f - 31n2) t + 0{t'^, 1/g) and fi{t) = 0{l/g). To determine the functions fo{t) 
and fi{t), we solve the quantization conditions (I2.10p and, then, use the obtained expressions 
for Cp{g) to resum the series in (12.80 . In this way, we obtain that fo(t) and fi(t) are given 
by a linear combination of the ratio of Euler gamma-functions (see Eq. ( ]A.4p in Appendix A). 
Their substitution into (12.151) yields the expression for functions T^lAngit) in the double scaling 
hmit 1261 



T^{4:ngit) + iT_{4:'Kgit) 



Vo^AT^gt)—— 



+ 



Ang 



n'i-t) 

T(i)r(i + t) 
Atra + t) Atra-t) 



r(|)r(i 



t) 



(2.16) 



+ 



Ang 



31n2 1 

4 ^ si 



mm 



t) 



t) 



r(|)r(i + t) 
stra + t) 



where ellipses denote sub leading terms suppressed by powers of 1/g and the functions Vo^Angt) 
and Vi{ATTgt) are defined in (I2.13p . The relation (I2.16P is consistent with the expansion (12 . 1 21) 
and it can be used to compute the 1/g correction to the mass gap (12. 5p . 
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We are now ready to compute the mass scale (I2.5p . To this end, we substitute (I2.16P into (12.50 
and work out the asymptotic expansion of the t— integral at large g. For the first term in the 
right-hand side of fl2.16l) this calculation was already performed in [7|. It leads to the expression 
for m which coincides with the mass gap (11. 4p found from the string theory considerations [6]. 
Taking into account the remaining terms in the right-hand side of (12.161) we should be able to 
compute the subleading correction to m. Calculation goes along the same lines as in [7j and it 
leads to (see Appendix A for more details) 



kg 



1/4 g-7rg 



1 + 



3 - 61n2 
327r^ 



+ 0{l/g' 



(2.17) 



with k = 2^/'*7r^/^/r(5/4). Comparing this relation with (12.11) we conclude that 



3 3 

mi = m 2 . 

32 16 



(2.18) 



This result is in an agreement with the numerical value found in the last reference in [25] . 

The coefficient in front of {3271 g)^^ in (12.171) is given by the sum of two terms of different 
transcendentality. We observe that ~ In 2 term can be absorbed into redefinition of the coupling 
constant 



9 



9 - 



31n2 

in 



(2.19) 



so that the mass gap (12.171) in terms of the shifted coupling g' looks as 



m 



i^9'[ 



r(5/4) 



1 + ^^ + 0(1/^ 
327r^' ^ 



(2.20) 



It is interesting to notice that similar simplification also occurs for the cusp anomalous dimension 
rcusp(5') at strong coupling. As was found in ^24], the expansion coefficients in the strong coupling 
expansion of Fcuspls') also involve terms ~ ln2 but they disappear after one re-expands the series 
in 1/g'. This suggests that the expansion parameter at strong coupling is g' rather than g. 



2.3 Induced renormalization scheme 



As was already mentioned, the mass scale m emerges in the AdS/CFT correspondence through 
dimensional transmutation mechanism in an effective two-dimensional theory describing dynam- 
ics of massless excitations in the AdSs x sigma model. As a result, the dependence of the 
mass scale of the coupling constant is dictated by the renormalization group. 

The coupling constant in the effective theory depends on the scale and is related to the 
coupling of the (conformal invariant) AdSs x sigma model as g'^ifJ^) = l/{2g) with the scale 
/i ~ 1 defined by masses of massive excitations [6l[28]. The coupling g{fi) satisfies the Gell-Man- 
Low equation 

dg 



and it leads to the following expression for a renormalization group invariant scale 



djjL 



A = jie 



rg dg 
J /3(S) 



2l3og^ 
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1 



1 + :^ ^ 



(^1 



/52 , . 



(2.21) 



(2.22) 
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Then, the relation ( I1.3P between the scahng function and 0(6) sigma model implies that j3{g) 
in (12.211) should coincide with the beta-function of this model. In bosonic two-dimensional 0(6) 
sigma model the beta-function coefficients are given by 



1 

TT 



Pi 



27r2 



(2.23) 



Notice that the beta function (12.211) depends on the renormalization scheme starting from 0{g'^) 
term. The same is true for the scale A while the mass scale m should be scheme independent. 
The two scales are related to each other as 



cA , 



(2.24) 



where the coupling independent factor c is needed to restore the scheme independence of m. In 
the special case of the MS scheme, this relation takes the form (ll.lSp with c^jg = e''. 

Replacing the beta-function coefficients in (I2.22p by their actual values (12.231) we obtain from 

(EMU 



m 



(2.25) 



Let us now compare this relation with the expression for the mass scale (12.171) obtained from 
exact solution of the FRS equation. We find that the two expressions indeed coincide upon 
identification 



1 



r(/i) = — 



r(5/4) 



P: 



(FRS) 



1 



(61n2 - 1) . 



(2.26) 



Here we introduced the subscript (FRS) to indicate that these expressions are valid in a particular 
renormalization scheme dictated by the FRS equation. 



3 Energy density in the two-dimensional 0(n) sigma model 

The aim of this section is to calculate the ground-state energy density of the 0{n) sigma model 
for large particle density p. The ground-state energy density eo(n) can be obtained either from 
the solution of the integral TBA equations, Eqs. (11.91) and (11.81) . or from perturbation theory. To 
compute the constant (12. 4p and, then, to make a comparison with the results of Refs. |28j and 
pO], one needs the second subleading correction to £o(6) in the large p limit analogous to ei in 
(12.31) . Since the techniques on the TBA side are not developed enough to perform the expansion 
at such depth, we determine this correction using standard perturbation theory. 

3.1 Perturbative calculation of the free energy density 

The fundamental fields of the 0(n) nonlinear sigma model are S{x) = {S^, . . . , S*") subject to the 
constraint Yli = 1. The theory has global 0{n) symmetry and the corresponding conserved 
charges can be written as 

Q'' = I Jl'dxi , = S'd^S^ - S^d^S' . (3.1) 
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We couple an external field to the conserved charge Q^^ and define the theory by its Euclidean 
two-dimensional Lagrangian: 

C{x) = ^ [d^Sd^S + 2th{S'doS^ - S^doS') + {{S^)^ + ■■■ + (5")^ - l} - 2u^S^] , (3.2) 

where A is the bare coupling constant and /;,— dependent terms are chosen in such a way that 
they modify the Hamiltonian of the model by term {—hQ^"^). 

In order to avoid infrared divergences we introduced in the right-hand side of (13.21) an extra 
term with regulator u which is going to be put to zero at the end of the calculation. This extra 
term fixes the classical ground-state to be 

= 1, = = ■■■ = S"" = 0. (3.3) 

We parameterize the small fluctuations around this ground-state by exploiting the remaining 
symmetries: 



= yi - A2(y2 + ;22)^ S' = Xy^ S' = Xz\ = A^"-^ (3.4) 

where the fields z = {z^, . . . , 2;""^) form the vector representation of the unbroken 0(n — 2). We 

substitute (13. 4p into (13. 2p and expand the Lagrangian to second order in the coupling A to get 

Cix) = A-2/:„2 + jCo + XCi + X^C2 + 0{X^) , (3.5) 
where the various terms depend on the parameter = h"^ + u"^ and look as follows 

= -^iM' + u;'), (3.6) 
£0 = ^ [d^yd^^y + d,,zd^z + + M^z^] , 



£1 = -ih{y' + z')d,y , 



1^2 = -{yd^y + zd^z){yd^y + zd^z) + —{y 



^2\2 



Here C-2 is just a constant, Cq defines the kinetic term for the two-dimensional fields y{x) and 
z{x) while Ci and C2 define, correspondingly, cubic and quartic interaction vertices. Notice that 
the Lorentz covariance of Ci is broken by the external field. 

Our goal is to calculate the change in the free energy density J^{h) — JF(0) relative to its value 
for /i = 0. The free energy density is defined as 



^-VT(h) 



j VyVz e-^'^°''^^''\ D = 2-e. (3.7) 



It is ultraviolet divergent and we used dimensional regularization to define it properly. Substi- 
tution of (13.51) into (13. 7p yields the perturbative expansion 



e-^^W = e"^^ / Pi/Pze-/'^"^^«(^) (3.^ 



X 



1-A^ 



J d^x£2(a;) + y J d^'xCiix) J d^x'£i(x') + O(A^) 
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Here the term linear in Ci is absent since it only involves odd powers of ?/— field and, therefore, 
vanishes upon integration. Taking the logarithm of both sides of (13.81) . dividing by volume 
V = J d^x and subtracting J^{Q) we get 

r{h) - ^(0) = {h) + {h) + A2^(i) {h) + O(A^) . (3.9) 

We calculate each term separately for a finite uj and put a; ^ at the end of the calculation. We 
relegate the details of the calculations to Appendix B and present here only the results. 

The first term in the right-hand side of (13. 9p describes the classical contribution to the free 
energy 

^(-^)(/.) = -^. (3.10) 

The next order term, J-'^^\h), sums up the quadratic fluctuations of y— and z— fields and it is 
related to the determinant of the kinetic operators in Cq. Using dimensional regularization it 
can be written as 

^(°)(/.) = ^/.-^ {7 + i + ' ^ = + ■ ^2-^^) 

Finally, going through calculation of the the third term in the right-hand side of (13. 9p we find 

At this point, we combine together three terms and obtain the following expression for the free 
energy density 

The first two terms in this expansion were already calculated in [12]. The result for the O(A^) 
term is new and it is needed to make the comparison with the results of Refs. [28] and |30j . 



3.2 Renormalization of the free energy density 

Ultraviolet divergences in the free energy (I3.13P can be removed in the standard way by expressing 
the free energy in term of the renormalized coupling and using the renormalization group to 
improve the result. In the MS scheme, the relation between bare coupling A and renormalized 
coupling g{fi) reads 

A^^(^eV^)'z,^, = + + (3.14) 

e e 

where 7 is defined in (13.111) and /5o,i,2 are the coefficients of the /3— function of the 0{n) model 
up to three loops 

/^^ = m = -Pof - Pif -f32f + ... (3.15) 
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In the MS scheme they have been determined in [31] to be 



^0-^' ^2 , (3.16) 

where we introduced the superscript to indicate that /?2 is scheme- dependent. 

Note that the couphng g introduced here is analogous to the couphng g introduced in Sect. 2.3. 
However the important difference between the two couphngs is that they are defined in two 
different schemes. Indeed, it is straightforward to verify that for n = 6 the beta-functions (12.231) 

and (I3.16P coincide up to two loops, but they differ starting from three loops, (3^^^^ ^ 
Still, the two couplings are related to each other through a finite renormalization. 

With the relation (I3.14p taken into account, the renormalized free energy density reads 

:Fih) - ^(0) = -y - 2/5o (in^ + - 2P^~g' (ln^ + + 0(^^)| . (3.17) 

We verify with a help of (I3.15P that the right-hand side of this relation does not depend on 
the renormalization scale /x. This suggests to express J^{h) — J^{0) in a renormalization group 
invariant way. 

It is important to keep in mind that, contrary to the free energy density, the coupling g{n) is 
not a physical quantity. We can explore this fact to define a new coupling constant to our best 
convenience. The running of the coupling is determined by the Gell-Mann-Low equation (13.150 
and it depends on the scale Ay^ given by a general expression ( 12.22^ in the MS scheme. Using 
this scale we define a new universal coupling a{h) as 

-+eina = ln-^, ^=^ = ^. (3.I8) 

Differentiating both sides of (13.181) with respect to h we find the coupling a defines the scheme 
in which beta-function is given to all loops by 

h-JJ- = = ■ (3.19) 
an 1 — 

It has the advantage that any other coupling constant depends only polynomially on a. In 
particular, for fi = h we find from (I3.15P and (13.180 



~9 (h) = ^^[a + ^a' + Oia'')j. (3.20) 

Then, perturbative expansion of the free energy density (I3.17P in a takes the form 

J'ih) - ^(0) = -Poh' - ^ - Y + 0(«')} . (3.21) 

Here the first two terms reproduce the result of pjj and the third one is the sought next order 
correction. 

Let us now determine perturbative expansion of the energy density £{p)- It can be obtained 
from the free energy density (I3.2ip through Legendre transformation 

e{p) = J^{h) - T{0) + ph , p = -j^'{h). (3.22) 
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Using (13.18P and (I3.2ip we find explicitly 



p^a{h) 
P = Woh 



a(h) ia^ih) ^. ^. 



(3.23) 



a{K) 



These relations define a parametric dependence of the energy density e on the particle density p. 
To express e entirely as a function of p we introduce yet another coupling a(p) defined by 



i + (e-l)lna = ln— ^ 



(3.24) 



Replacing p by its expression (13.231) and making use of (13.181) . we establish the relation between 
the two couplings 



&{p) 



+ (e-l)ln«(p) 



a{h) 



+ {i-l)\Yia{h) + 0{a^ 



(3.25) 



leading to a{h) = a{p) + O(a^). 

As a result, we obtain from (I3.23P the energy density in the 0(n) model as 

e(p)= Ae|a + ^ + ^ + 0(a^) 



(3.26) 



According to (I3.24p . the coupling a depends on the scale A^. To make a comparison with the 
string theory calculation we need its expression in terms of the physical mass m. To this end, we 
take into account the known relation between the two scales 



C = ln 



m 



A 



(31n2-l)e-lnr(l + 



(3.27) 



MS 



and express the coupling (I3.24p in terms of p/m as 

- + (^-l)lna = ln^ + A 

a m 



A = C + \n (vrO 



(3.28) 



with ^ = l/(n — 2). Being combined together, the relations (13.260 and (I3.28P define the first few 
terms of the perturbative expansion of the energy density in the 0(n) model in the perturbative 
regime p 3> m. 



3.3 Numerical analysis of the TBA equation 

As shown in the second paper in [12j, the same function 6{p) can be obtained from the solution 
of the integral TBA equation (II. 9p in a parametric form given in (II. 8p . Using the generalized 
Wiener-Hopf technique one can transform equation (II. 9p to a form more suitable both for finding 
the large p asymptotic expansion of the energy density and for numerical calculation. Employing 
this technique, Hasenfratz, Maggiore and Niedermayer [12] calculated the first two terms in the 
large p expansion of e{p) and computed the ratio m/A^jg. 
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Since for the calculation of the next order term we have to leave the beaten paths we decided 
to check the result of the perturbation theory numerically for the 0(6) model. The precision mea- 
surement would require to perform several integral transformation, but a rough estimate can be 
obtained directly form the original TBA equations (11.81) and (11.91) as follows: First the parameter 
B is fixed as an integer in the range B = 5, ... ,14. On the interval [—B, B] the pseudo-energy 
x{0) is discretized on 2^^ -^ 2^^ points. Then the integral equation ( 11. 9p is solved by iteration. 
(Technically, we used the fast Fourier transformation routine of the numerical programming lan- 
guage octave to perform the convolution). From the solution, x(6', 5), we calculated the ratios 
p{B)/m and e{B)/w? with a help of (II. 8p . Our relative precision for these quantities was as 
good as 10~^. Then we calculated the coupling a{p) by numerically solving the equation (I3.28p . 
The results are displayed in Table [H 



B 


5 


6 


7 


8 


9 


a 


0.23355(3) 


0.19097(5) 


0.16135(3) 


0.13958(8) 


0.12293(8) 


p/m 


94.8911(7) 


286.651(8) 


850.661(7) 


2492.50(2) 


7233.93(2) 


e/m^ 


1856.75(4) 


13560.90(5) 


99410.4(0) 


730337.7(8) 


5373044.(2) 


B 


10 


11 


12 


13 


14 


a 


0.10979(9) 


0.099173(3) 


0.090405(5) 


0.083049(9) 


0.076792(3) 


p/m 


2.08399(8) ■ 10^ 


5.96823(1) ■ 10^ 


1.70101(3) ■ 10^ 


4.82856(7) ■ 10^ 


1.36601(2) ■ 10^ 




3.95663(6) • 10^ 


2.91555(6) ■ 10® 


2.14945(9) ■ 10^ 


1.58524(9) ■ 10^° 


1.16947(3) ■ 10" 



Table 1: Numerical result for the energy and particle density in the 0(6) sigma model as function 
of parameter B. Unreliable digits are displayed in parenthesis. 



Finally we fitted the ratio eji^fp'a) to the expression 

e Tx ( a 0? ^ , oA 

?5i = l(i + I + T + «'*')' 

which follows from f l3.26p in the case of the 0(6) model. Feeding in the known first two coefficients 
we obtained 0.12(4) for the coefficient of the interesting O(a^) term in the right-hand side of 
f l3.29p . This is in good agreement with the predicted value |, recalling that here we are dealing 
with the coefficient of a second subleading correction. 

To visuahze the result we collect the data on Figure 1, where we plot e/^p^a) against the 
coupling a. Clearly the sum of all the three calculated perturbative contributions approaches 
nicely the numerical solution giving a convincing confirmation of both. We obtained similar 
result for the 0(3) model as well, thus confirming our perturbative calculations numerically. 

3.4 Large j expansion of the energy density 

To make a decisive comparison with the results of Refs. [2B] and [20], we have to examine the 
behavior of the energy density e{p) for large particle density p ^ m, or equivalently j ^ m with 
the scaling variable j = 2p defined in ( II. 3p . 
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Figure 1: e/{p'^a) is plotted against a. Boxes represent the numerical TBA data with invisible 
numerical errors. The lines correspond to perturbative corrections as marked on the figure. 



In this region, the coupling constant fl3.28p is small and it has the form 

1 7 7 In In — ua 

- = In — + ai Inln — + as + as + t- + ... . 

am m In — In — 



The coefficients a^ can be determined recursively from (13.281) leading to 

ai=^-l, a2 = A, a3 = (^-l)^ a4 = (^-l)A. 



(3.30) 



(3.31) 



with the constants ^ and A defined in (I3.18P and (13.281) . respectively. We substitute the relation 
(I3.30p into (13.261) and obtain, after some algebra, the following expression for the energy density 
of the Oin) model 



^O(n) 



J ^ (1-0 
In^ 



41n^ 

m 



, (1-0^ 

In^i 

m 

where the notation was introduced for 



In ( K In — I H — 
' m 2 



+ . . 



In K 



A _ (|-31n2)e + ln(r(0/vr) 



1-e 



1-e 



(3.32) 



(3.33) 
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We recall that the parameter j = 2p defines the density of particles with mass m and the relation 
fl3.32p only holds for j m. 

We are now ready to perform a comparison of the energy density of the 0(6) model and the 
scaling function (12 .31) in the AdS/CFT. For n = 6 the constants fl3.33p take the following values 

^ = \, ln« = l-iln2-^lnr(|). (3.34) 

Taken into account these relations we find that, in agreement with (11. 3p . the expression for 2£:o(6) 
indeed coincides with the two-loop result for the scaling function (I2.3P with 



5(1 

^1 — "7 



7 

n=6 36 



(3.35) 



(1-0^ 

Finally, we substitute the relations (12.180 and (13.350 into (12.40 and compute the constant go2 as 

q^^ = ^ + ^-^l,,2 + --- = -~-\n2. (3.36) 
^" 8 4 2 2 36 4 2 ^ ' 

Comparing this relation with (II. 7p we conclude that our result for go2 is in agreement with the 
result obtained from the quantum string Bethe ansatz [30] . 



4 Conclusions 

In this paper, we applied methods of integrable models previously developed for the two-dimen- 
sional 0(n) sigma- model to study the scaling function in the AdS/CFT. The starting point of 
our analysis were the relations (II. 3p and (II. 4p which relate the energy density of the 0(6) model 
to the scaling function and identify the mass gap of the 0(6) model m with a new dynamical 
scale in the AdS/CFT. 

These relations are extremely nontrivial given the fact that the scale m has a different origin 
in the models under consideration. The 0(6) model is asymptotically free at short distances and 
the mass scale arises due to a nontrivial dynamics in the infrared. At the same time, Af = 4 SYM 
and string sigma-model on AdSs x S^ are conformal invariant theories and, therefore, they do 
not generate any scale. In the string theory, the scaling function describes quantum corrections 
to the folded string spinning on AdSs x S^ and it is the underlying classical configuration that 
introduces the mass scales. In a similar manner, in A/" = 4 SYM theory the scaling function 
follows from the analysis of the Bethe ansatz equations in the limit (II. ip . In this limit, the Bethe 
roots condense on the real axis and their distribution depends on the parameters (11.10 fixed by 
the quantum numbers of Wilson operators. It is therefore remarkable that our calculation of 
the mass scale in A/" = 4 SYM theory reproduces the known result for the mass gap in the 0(6) 
model in the special renormalization scheme dictated by the FRS equation. 

For g ^ oo and j/m = fixed, the scaling function can be found exactly by solving the TBA 
equations for the energy density of the 0(6) model. Finding solution to these equations in the 
perturbative regime j ^ m beyond the leading order proves to be a difficult task. We demon- 
strated that the problem can be circumvented by employing methods of standard perturbation 
theory to calculating the free energy of the 0(n) model in the presence of an external field. In 
this way, we computed the two-loop correction to the energy density of the 0(n) model and 
verified that it agrees with the numerical solution to the TBA equation. Then, we combined 
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this correction (for n = 6) with the two-loop correction to the mass scale m and calculated the 
scaling function. 

Comparing the obtained expression with two different predictions for the scaling function 
coming from two-loop stringy corrections to the folded spinning string solution [2^ and from 
thermodynamical limit of quantum string Bethe ansatz equations |lI51 |29l [30], we found an 
agreement with the latter one. The agreement should not be surprising since the FRS equation 
was originally obtained from the all-loop Bethe ansatz for the dilatation operator in = 4 SYM 
theory in an appropriate scaling limit. 

The question remains however what is the reason for a disagreement between our result for 
the scaling function and the explicit two-loop stringy calculation. Notice that the difference 
only amounts to a constant term (see Eqs. ( II. 6p and 01. 7p ) while logarithmically enhanced terms 
coincide. This implies that the two-loop stringy result is consistent with the relation (11.31) between 
the scaling function and the energy density of the 0(6) model but it leads to the expression for 
the two- loop correction to the mass scale (12.11) which is different from (12.181) . 

mf = -— - — ln2- ^K. (4.1) 
^ 32 16 4 ^ ' 

This suggests that either the all-loop Bethe ansatz does not predict correctly the mass scale, or 
the two-loop stringy calculation needs to be revisited. At present stage we can not discriminate 
between these two scenarios and the question requires further investigation, [f] 
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Appendix A: Calculation of the mass scale 

Let us now compute the mass gap ( 12. Sp . At large g the integral in ( 12. 5p receives a dominant 
contribution from t ^ g. In order to evaluate ( 12. 5p it is convenient to change the integration 
variable as t — > Airgit. Then, we get from ( 12. 5p 

8\/2 

m=^e-^f+Am, (A.l) 

potential difficulty in comparing the two predictions is that the string computation |28j was performed in a 
scheme in which two-loop beta-function coefhcient is zero. This scheme is related to the MS scheme by a singular 
coupling redefinition and it was previously used [32] in the studies of two-dimensional bosonic 0(n) model. We 
are grateful to A. Tseytlin for drawing our attention to this fact. 
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where Am is given by integral of r{4TTgit) along the imaginary axis 



Am 



TT 



dt e'4^9*-*^/4 



t 



{T^{Angit) + iV _{A7!-git)) + c.c. 



According to (12.161) the functions r±{47T git) have the following form 

T^{4ngit) + iT_{4ngit) = Vo{47Tgt)fo{t) + Vi{4ngt)fi{t) 



(A.2) 



(A.3) 



where /o and /i are defined as coefficients in front of Vq and Vi, respectively, in the right-hand 
side of (12.161) . They are given by a linear combination of the ratio of Euler gamma-functions 



m 
flit) 



r(|)r(i 



Ang 



'3 

r(i)r(i 



A-Kg 



31n2 1 

4 ^ 8t 



r(|)r(i 



t) r(i)r(i + t) 



+ 0{g 



-2\ 



r(|)r(i 



t) 



Atva + t) 



Atva-t) 



+ 0{g- 



(A.4) 



Notice that fi{t) is suppressed by factor l/^A-ng) compared to foif). 

To work out the large g expansion of the integral ( 1A.2|) it is convenient to use the Mellin- 
Barnes representation for the functions Vo(47r(7t) and Vi{ATigt). Using the integral representation 
(Km we find 

^ ' -j){Angty I du{l-uy-^/\l + uy+^/^ (A.5) 

' —5—ioo 



Vo{Angt)e- 



TT 



2TTi 



-1 



and similar representation also exists for Vi{A7cgt) e ^^^3*. Their substitution into (1A.3I) and ( ]A.2I) 
yields after m— integration 



Am 



^V2g 



r(-j)r(j + 1) 



—5—ioo 



2ni 

dt {Aty 



r(j + 2) 



(27r^)^e-*"/^ 



(A.6) 



-fo{t) + hmj + l) 



+ c.c. 



To find the asymptotic expansion at large g we deform the integration contour in the complex 
j plane to the left and pick up the contribution of poles at negative j. These poles come from 
r(j -|- |) and from t— integral. Let us start with the latter ones. 

By definition, fo{t) and fi{t) are real meromorphic functions of t with simple poles located at 
t = ±1, ±2, .... Then, integration over small t produces simple poles located at negative integer 
j. However, due to the presence of l/r(j + 2), all these poles except j = —1 produce vanishing 
contribution to the j— integral. Calculating the residue at j = — 1 we find 



Am = e-^ T{m-ly2ng)-' e^ i/o(0) 



8^2 



,„ , „, +c.c. + ... = --^e--^ + ..., (A.7) 

where ellipses denote the contribution of poles produced by r(j ' + |) in flA.6l) . We notice that 
( ]A.7p cancels against similar term in the right-hand side of (lA.ip and, therefore, m is determined 
by the contribution of poles at j = — |, — |, . . .. Going through calculation of residues we obtain 



m 



Aig 



71^ 



-3/4 g-7rg 



—5—ioo 



foit) + Iflit) - i^ngty'^oit) + Oig-') 
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where the integration contour goes to the left from the branch cut that starts at t = 0. Calculating 
this integral and taking into account analytical properties of the functions (1A.4I) we find 



m 



1 3 

/o(-i) + 2/i(-i) + ^ 



fo{-l) + Oig 



-2\ 



leading to 



21/2 



61n2 - 3 
327r^ 



0(9 



-2a 



(A.9) 



(A.IO) 



in an agreement with fl2.17p . 

The strong coupling expansion of the mass scale (lA.lOp can be systematically improved by 
taking into account subleading 1/g corrections to the functions /o(t) and in (]A.4|) . Assuming 
that higher order corrections do not modify analytical properties of these functions and taking 
into account a contribution to (1A.6I) from an infinite sequence of poles at j = —3/4 — n (with n 
nonnegative integer) we find after some algebra 



m 



(A.ll) 



Here the notation was introduced for the functions 



Uoiy) 



1 + 



32y 



+ . 



'3^ 
2^ V4y 



dt e-^yH~'/\i + tY/^ = (2|/)-3/4r(|) 

U-{y) = - dt e-^yH-^'\l + t)-3/4 = {2y)-^l%T 
2 JO 

which can be expressed in terms of Whittaker functions of second kind 

U^{y) = \T{^y-^eyWy,,„2{2y), 
U^{y) = \T{W^y)-"'eyW^y,,,{2y). 



32y 



+ 



(A.12) 



(A.13) 



Appendix B: Perturbative calculation of the free energy 
density 

In this appendix we calculate the constant {0{\^)) and quadratic {0{\^)) terms in the pertur- 
bative expansion of the free energy (13.91) . 

The constant term sums up the quadratic fluctuations of fields y{x) and z and it is given by 
a logarithm of the ratio of determinants of the kinetic operators: 

^{o)(/,) = "^-Zl [Trlog(-a2 + M^) - Trlog(-92 + cu2)] , (B.l) 

with = h? + uj"^ and V being the volume factor. Fortunately computing the difference also 
provides a regularization. By differentiating and integrating with respect to M and uj we can 
cast the result into the form 

^(o)(/,) = liin r dm'"" I Z"",^, , ^ , . (B.2) 

^ ' 2 <.^-*oJ^2 J (27r)^p2 + ^/2 ^ ^ 
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The momentum integration can be performed in the dimensional regularization with D = 2 — e 
as 

-+7 + 0(e) 
e 



i m = / tt:^; t: = m — 



(B.3) 



with 7 = r'(l) + ln(47r). As we will see in a moment, the same integral appears in higher order 
calculations. Using fIB.Sp . we find for the constant term (1B.2P (up to corrections vanishing as 

6-0) 

^(0)(/.) = !^/.-4l + 2 + a. (B.4) 



47r l^e 2 2^ 

The 0{\^) term in the expansion fl3.9p describes the two-loop correction to the free energy. It 
receives contribution from the last two terms inside the square brackets in (13. 8p 

T^'\h) = j d""^ (A(^))o -^jd^'^j d^'x' {C,{x)C^{x')), , (B.5) 

with Ci{x) and £2(2;) defined in (13.61) . Here the subscript '(0)' indicates that the expectation 
values are evaluated with the measure / VyVzexp {-J d^xCo{x)) (see Eqs. ([31]) and ^M ). 
The expression for f d^x {C2{x))o = {€-2)0 can be obtained in terms of the two VEV's 

{y\ = I{uj) , {z\ = {n- 2)/(M) . (B.6) 



Using translational invariance (9^(?/^))o = and the equation of motion d'^y = uj'^y, we get 
{{yd^yY)o = —^{y^)o = —uj'^{y'^)1. Making use of similar relations for z— field together with the 
factorization property {yd^yzd^z)Q = {yd^y)Q{zd^z)Q = 0, the contribution of £2 to (IB.Sp can 
be written in the following way: 

(A)o = -p'H-^-^^^lHM) (B.7) 

+ — [SPico) + 2{n - 2)I{lu)I{M) + 2{n - 2)P{M) + {n - 2fP{M)\ . 
8 

We remove the infrared cut-off by taking the — limit and obtain the relevant contribution 
from £2 as 

1 /" n —2 

- j rf^x (£2(x))o = (£2)0 = i^h^l\h) . (B.8) 

The contribution to the free energy ( 13. 8p quadratic in £i(x) has the form 

^jd^'xj d^x'(£i(a:)£i(a:'))o = \ j d''^ (£i(a;)£i(0))o . (B.9) 

By writing 

£1 = -ih{y^ + z^)d^y = -ih {\d^y^ + z^doy) , (B.IO) 

we can see that the first term, being a total derivative, can be dropped. The remaining term 
gives rise to the diagrams shown on Fig. [2l The contribution of the first diagram on Fig. [2] 
involves the factor ~ J d^x {doy{x)doy{0))o with the integrand being a total derivative again. 
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Figure 2: Two- loop diagrams contributing to (|B.9p . Solid lines denote two-dimensional scalar propa- 
gators. 



Thus the right-hand side of flB.9p only receives contribution from the second diagram on Fig [2j 
Taking into account all possible contractions one finds 



-Jd^x{C,{x)Cm)o = -h\n-2)Jd^Xj ^^^^n^^ + M^J (2.)^ p^ + J {2n)^ + ^ 

(B.ll) 

where the factor [n — 2) counts the number of 2— fields circulating inside the loop. Doing the 
p— integration we observe that despite the fact that the integrand is not Lorentz covariant, the 
integral in the right-hand side of ( IB. lip could only depend on and uj"^ and, therefore, it should 
be Lorentz invariant. This allows us to simplify the p— integral as: 



^2 „ipx 



d p Pq e 
(27r)«p2 + ^2 



D J {2n)^p^ + uj^ D ^ ' ^ ' 



(B.12) 



As a consequence, the a;— integral in flB.llP becomes trivial and the remaining momentum inte- 
gration gives P{M). Then, taking the uo ^ Q limit, we find from (IB. IIP 



^ j rf^x (£i(x)/:i(0))o = -^(n - 2)l\h) . 



(B.13) 



Substituting the relations (IB.SP and (IB.lSp into (IB.SP we derive the sought second order correction 
to the free energy density 



T^^\h) = h''{n-2)I^{h) 



1 

'd 



n 



2^2-2. ^ 1 



+ 7 + 



(B.14) 
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